Transient fluctuation relations for time-dependent particle transport 
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We consider particle transport under the influence of time- varying driving forces, where fluctuation 
relations connect the statistics of pairs of time reversed evolutions of physical observables. In 
many "mesoscopic" transport processes, the effective many-particle dynamics is dominantly classical, 
while the microscopic rates governing particle motion are of quantum-mechanical origin. We here 
employ the stochastic path integral approach as an optimal tool to probe the fluctuation statistics 
in such applications. Describing the classical limit of the Keldysh quantum nonequilibrium field 
theory, the stochastic path integral encapsulates the quantum origin of microscopic particle exchange 
rates. Dynamically, it is equivalent to a transport master equation which is a formalism general 
enough to describe many applications of practical interest. We apply the stochastic path integral 
to derive general functional fluctuation relations for current flow induced by time-varying forces. 
We show that the successive measurement processes implied by this setup do not put the derivation 
of quantum fluctuation relations in jeopardy. While in many cases the fluctuation relation for a 
full time-dependent current proflle may contain excessive information, we formulate a number of 
reduced relations, and demonstrate their application to mesoscopic transport. Examples include 
the distribution of transmitted charge, where we show that the derivation of a fluctuation relation 
requires the combined monitoring of the statistics of charge and work. 
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PACS numbers: 05.40.-a, 05.60.-k, 73.23.-b, 72.10.-d 



I. INTRODUCTION 

Fluetuation relations (FRs) have recently emerged 
as a new powerful set of concepts in statistical 
physicsji"— They relate the stochastic fluctuations of 
systems far from equilibrium to their dissipative prop- 
erties, thereby generalizing the well-known fluctuation- 
dissipation theorem;^ and provide ways to quantify 
the degree of irreversibility of nonequilibrium processes. 
The possibility to formulate exact statements (some- 
times referred to as "fluctuation theorems") about 
generic nonequilibrium systems - both classical2i2ii2riii 
and quantum^ d^d'^ - has ignited a burgeoning research 
activity. On the experimental side, first tests have al- 
ready appeared, e.g., for soft matter—"— or mesoscopic 
systems ii2r— 

The concept of FRs is frequently applied to the statis- 
tics of variables of thermodynamic significance, e.g., 
work, heat, or entropy. This is exemplified by the Crooks 
rclatioE^iSiii 
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where P{W) is the probability that an amount of work 
W is done on a system in a given driving protocol, i.e., 
when an external time-dependent force ft is applied to 
the system during a time interval t e [— t, r]. Accord- 
ing to the Crooks relation, the ratio to the probability 
Pb{~W) of negative work done during the "backward" 
protocol, i.e., when the time- inverted force /_t acts on 
the system, is given by a Boltzmann-type factor, where 
/3 — T~-^ is the inverse temperature (fcs = 1 throughout) 
associated with the initial equilibrium state, and AF is 



the thermodynamic free energy difference between final 
and initial state. Relations of this type may be applied 
to gain access to thermodynamic data (e.g., AF) from 
fluctuation statistics. From Eq. (jl.l|) one obtains the 
celebrated Jarzynski equality;^ 
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where the average is over all process realizations under 
the force protocol ft. On average W > AF, which means 
that the "sum rule" (|1.2|) controls the cumulative weight 
of rare events. Relations of this type provide rigorous 
bounds on the behavior of thermodynamic observables. 
Fluctuations around these bounds have been analyzed 
recentlyj^"— 

The application of external forces in many-particle 
systems generally leads to transport, with current flow 
and re-adjustment of particle concentrations. For time- 
dependent driving forces, FRs for the ensuing current 
profiles carry particularly rich information. Here, the 
notion of "transport" is to be interpreted in a very gen- 
eral sense: it may refer to the changing number of indi- 
viduals in a biological quasi-species model, to the elec- 
tric current flow in a mesoscopic conductor, to the num- 
ber of agents in a chemical reaction, etc. Unlike with 
FRs for "global" (i.e., integrated over time) variables, 
the time-resolved information on transient current flow 
is stored in a time-dependent function / = {It}, and 
one needs to study functional probability distributions 
P[I], rather than functions like P{W). To be specific, 
transport through a "system" exchanging particles with 
M "reservoirs" is described by currents Iu=i,...,m flowing 
out of the system, and P[I] = P[Ii, . . . ,Im]- Impor- 
tantly, due to the discreteness of particle exchange with 



the reservoirs, the current flow enhances the noise level of 
the system. Far off thermal equilibrium, this "shot noise" 
often becomes the dominant source of fluctuations, and 
the self-consistent description of the feedback cycle of 
currents generating noise and noise affecting the current 
flow becomes an important issue. 

Given the present interest in time-dependent meso- 
scopic transport phenomena, see for instance Refs. \24 - 
\2di . or in the work statistics under a quantum quench;^ 
it is important to extend the general constraints imposed 
by transient FRs to the quantum setting. The above ar- 
gument, however, needs to be applied with care in the 
quantum case. In fact, the existing literature on FRs 
appears to be essentially divided into a classical and a 
quantum part. This concerns both the systems under 
consideration as well as methodological aspects. Specifi- 
cally, quantum theories of FRs^ '^^i^'^'^^" — mostly rely on 
quantum measurement processes performed at the begin- 
ning and at the end of the protocol. This prescription is 
not directly suitable to transient situations, where a con- 
tinuous readout of, say, currents is required. The lack 
of commutativity of current operators at different times 
then becomes an issue, and the construction of a general 
quantum theory of transient FRs may seem a difficult 
task. Perhaps surprisingly, the derivation of a FR for 
quantum current flow goes through in unaltered form as 
long as the external driving forces vary on time scales cor- 
responding to classical frequencies. We will discuss this 
point in some detail below (see also Ref . |31I) . The situa- 
tion becomes particularly transparent in the many cases 
where the quantum system of interest actually operates 
close to the semiclassical limit: involving the dynamics 
of many particles, the action scales relevant to transport 
arc usually much bigger than h. In this case, the quan- 
tum Keldysh functional^"— stays close to its classical 
limit, the stochastic path integral^"— and the formal 
lack of commutativity of current operators ceases to be 
an issue. 

In an important early work,— Bochkov and Kuzovlev 
(BK) have formulated a general classical FR for current 
flow. BK relied on a symmetry analysis of the Markovian 
operators generating the dynamics of the system (see 
also the discussion in Ref. |30)- In contrast, our deriva- 
tions below are based on a path integral representation, 
which permits to explore the nonequilibrium fluctuation 
statistics of observables beyond the rigorous bounds im- 
posed by FRs. Also, it stands to reason that the BK 
functional relation is too general to be useful in applica- 
tions. However, the general result can be used to obtain 
more manageable derived fluctuation relations. For ex- 
ample, rather than probing the full profile / = {It}, one 
may consider the total charge transmitted into the i^th 
reservoir, Q^ = J_ dt I,j^f This is arguably one of the 
most important global variables characterizing a trans- 
port process. Under stationary transport conditions, the 
forward and backward protocols coincide, Pb{Q) = P{Q), 
and a Crooks FR for charge has been stated in the con- 



text of mesoscopic transport 
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This relation imposes nontrivial conditions on the gen- 
erating function for the full counting statistics (FCS) 
of charge transport4^ In particular, it implies that 
current cumulants of different order must be linked 
togetheri ^^i^^'^^" — Equation p.3p has also been probed 
expcrimentall y^^i^° in mesoscopic circuits using a quan- 
tum point contact charge detector4^ While the case of 
unidirectional single-electron counting is directly accessi- 
ble to experiments, recent progress has also been reported 
for bidirectional count ing<^2ii^ All these results apply to 
stationary regimes. Below we will show that in the exper- 
imentally relevant case of time dependent forces ft dras- 
tic things happen: the charge FR (|1.3p actually breaks 
down, but a more general FR for the joint probability 
P{Q, W) can still be formulated, see Eq. (|5.2p below. 
This "generalized Crooks relation" also ensures the va- 
lidity of Eq. (|l.ip for the statistics of work alone. Further- 
more, we will use the BK relation to derive cross-relations 
between nonlinear AC response coefficients, which can be 
put to an experimental test. 

The structure of the remainder of this article is as 
follows. In Sec. |TT] we discuss classical nonequilibrium 
transport processes in terms of master equations and 
the stochastic path integral. We derive the functional 
FR for currents under a transient driving protocol in 
Sec. mil The connection to the quantum theory is stud- 
ied in Sec. IIVI In Sec. |V] we discuss derived fluctuation 
relations and compare them to numerical simulations for 
a mesoscopic RC circuit. Some concluding remarks can 
be found in Sec. IVIl Various details of our calculations 
have been relegated to several Appendices. 



II. STOCHASTIC PATH INTEGRAL 

A. Master equation 

We are interested in the statistical properties of parti- 
cle currents flowing through a system in contact with M 
reservoirs. The probability Pt{n) that the system con- 
tains n particles at time t is given by the convolution 



Pt{n) = ^ Pt{n\n^r)p{n-T), 



(2.1) 



where Pt{n\n-r) is the conditional probability to evolve 
from an initial state n_7- at — r to n at t, and the 
weight p(n_T-) describes the probability of the initial 
state. Without much loss of generality, we take p to be 
of Boltzmann form, 



p(n)=e-''(^(")-^). 



(2.2) 



where U{n) determines the internal energy of the system 
and F — — Tln^^ exp(— /3t/(n)) is a free energy. The 



Markovian time evolution of Pt is governed by a one-step 
master equation;^ 



dtPtin) = -HgPt{n), 

M 



(2.3) 



H,{n,P) = ^^(l-eT^)g±,(n), 



=1 ± 



where the exphcitly time-dependent rates g+ {%) con- 
trol the flux into (out of) the system. The operator e^ 
(e~^) raises (lowers) n by one unit, i.e., we have the com- 
mutator \p,n] = 1. We assume that the rates obey the 
detailed balance condition 



9v,t 



K,^t{n)=dnU{n)~U,t, (2.4) 



where the functions f^, describe external driving forces 
on the system. For concreteness, we consider a cyclic 
protocol: starting from an equilibrium situation at time 
t = —T, fv,~~T = 0; the time-dependent protocol f^.t 7^ 
eventually ends at fy t = 0. The generalization of the for- 



J 



malism below to several types of particles, or to multiple- 
step master equations is straightforward. However, the 
generalization to situations where the reservoirs are at 
different temperatures requires a more substantial exten- 
sion of the formalism. In this case, one has to account 
for the energy transfer necessarily accompanying parti- 
cle transfer. We will briefly comment on this point in 
Sec. UVCl 



The condition (|2.4p is less restrictive than it might 
seem at first sight: it states that the logarithmic ratio 
of rates is governed by a cost function, T\i\{g^ / g~) = 
-~[Eu{n -|- 1) — Ey{n)\, measuring the difference in "ener- 
gies" Ei,{n) = —nfu + U{n) before and after a particle has 
entered the system through terminal v. Note that Ey {n) 
contains contributions linear in the driving parameters 
and the particle number, i.e., the driving couples to the 
energy balance of individual particles and not to particle 
interactions. With U{n) introduced in Eq. (|2.2p and the 
notation 9„C/(n) = U{n + 1)- U{n), we obtain Eq. ([2^ . 
Particle interactions then correspond to nonlinearities in 
U{n). 



System 


Variable n 


U{n) 


U,i 


electronic circuitsiii 


charge 


charging energy 


bias voltages 


molecular motors^ 


mechanochemical state of motor protein 


load potential 


ATP concentration 


chemical reaction networks^ 


number of reaction partners 


internal energy 


chemostat concentrations 


adaptive evolution^ 


allele frequencies 


log equilibrium dist. 


fitness gradients 



TABLE I: Examples of application fields for the master equa- 
tion. 
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Table U lists several application fields of present inter- 
est where the above model of discrete transport applies 
with little or no modification. For later reference, let us 
introduce one of the examples above in some more de- 
tail: consider a mesoscopic RC circuit, where the system 
corresponds to a central node ( "quantum dot" ) with n 
electrons held by it. The role of the reservoirs is taken 
by M = 2 voltage sources connected to the dot through 
resistors Ri and i?2- The driving of the system by a 
time-varying bias voltage Vt and its internal energy are 
given by /^.t = {-Y+'^eVt/2 - eV^s (with :/ = 1, 2) and 
U{n) = [n — l/2)-^Ec, respectively. Here, eT4ff is the 
effective chemical potential on the dot that needs to be 
determined self-consistently, Ec = e^/(2C) is the capac- 
itive charging energy of the dot, and we have taken an 
offset charge corresponding to Coulomb blockade peak 
conditionsi^ (In the rest of the paper we will set e = 1.) 
In earlier studies,— i^"— circuits of this type have been 
discussed within the framework of Langevin equations, 
where the dominant source of fluctuations was thermal 
noise. In contrast, we wish to include the more general 
mechanism of noise self-generated by transport. This 
physics can be described by the master equation (|2.3|) 



with sequential tunneling rateS)^^ 

1 ±Ky^tin) 



g^.M) 



i?^ e±/3«-.*(") -1' 



(2.5) 



with Kij,t in Eq. ()2.4|) . The Bose-Einstein function in 
Eq. (|2.5|) indicates a degree of quantum-mechanical in- 
put, to be discussed in more detail in Sec. IIVCI The 
rates (|2.5p comply with the balance relation (|2.4p . In ad- 
dition, it is straightforward to perform numerical simula- 
tions for the dynamics generated by the master equation, 
see Sec. El 



B. Path integral and current statistics 

To obtain the stochastic path-integral 
representation^SiiS of the above master equation, 
we allude to quantum-mechanical analogies and notice 
the similarity of Eq. (|2.3p to an (imaginary-time) 
Schrodinger equation. The unit-normalized quantity 



Z 



EPrir 



(2.6) 



can thus be written as imaginary-time path integral;^ 
Z= /7:>(n,p)e^-''*(P^*"-^«("'P»p(^-r). (2.7) 



The integration in Eq. (j2.7p is over smooth paths (n,p) = 
{(nt,pt)}, where the auxiUary "momentum" "pt G zR is 
integrated over the imaginary axis and the information 
on the actual discreteness of the evolution of nt is en- 
coded in the "Hamiltonian" Hg. 

In order to extract information on current flow, we 
need source fields. Specifically, time-dependent counting 
fields Xi^.t probing the current flow ly^t can be introduced 
by generalization of the Hamiltonian;^ 

i7,(n,p) ^ Fg(n,p,x) ^ E (l - eT(P-'x.)) ,,±(„). 

v.-k 

(2.8) 
Cumulants of the currents can then be obtained by func- 
tional differentiation, 

%Q to 



^Z[x] = 



^XutM 5xi^2M 



-iY:^!l.dtXuI. 



x=o 



(2.9) 



where Z[x] is the generating functional, see Eq. (|2.1ip 
below. Heuristically, the identification ^ ^ I follows 
from the fact that Xv enters the theory precisely like a 
vector potential. Much like in quantum mechanics, dif- 
ferentiation with respect to (w.r.t.) these vector poten- 
tials generates currents. More rigorously, the connection 
follows from a "Ward identity" of the functional inte- 
gral: temporarily considering the case of identical count- 
ing fields, Xv = X: the above functional derivative ob- 
tains j^\nZ ~ "^^{Iv^t), which we tentatively identify 
as the total current out of the system. On the other 
hand, x may be gauged out of the Hamiltonian by a shift 
p -^ p + ix, 8tt the expense of an extra term i J dt x^tn in 
the action. Differentiation w.r.t. x in the shifted repre- 
sentation obtains j2_in.Z[x] = —{dtUt). The equality of 
the two representations yields the continuity equation. 



(a,n)+^(/„)=0. 



(2.10) 



This shows that the differentiation w.r.t. x obtains the to- 
tal current, while differentation w.r.t. the fields Xi' yields 
currents through individual interfaces. 
The functional partition sum is 

Z[x] = j D{n,p)e-^o^^'P'^^p{n-r), (2.11) 
Sg[n,p,x] = -/ dt{pdtTi- Hg{n,p,x)), 



with the Hamiltonian (|2.8p . According to standard rules 
of probability theory, the probability distribution of cur- 
rents follows by functional integration over all x. 



P[I]= j Dxe'^^i-^^'^-'-Zix] 



(2.12) 



The discreteness of particle transport is encoded in 
the exponential dependence of the Hamiltonian on the 
"phase space momentum" p. Keeping this information 
is crucial, e.g., to properly resolve the statistics of rare 
events. When the discreteness does not play an impor- 
tant role, an expansion of e^ to quadratic order in p 
may be justified. This reduces the stochastic path inte- 
gral to the Martin-Siggia-Rose functionali^"— One may 
then continue to either integrate over momenta, which 
leads to the Onsager-Machlup path integral , ^^'^^ or de- 
couple the quadratic momentum dependence by an aux- 
iliary "noise field" , which generates an effective Langevin 
descriptioui^^i^ In the next section, we will employ the 
stochastic path integral p. lip to (re) derive a number of 
general FRs. 



III. TIME REVERSAL AND TRANSIENT 
FLUCTUATION RELATION FOR CURRENTS 

We proceed to derive a variant of the BK fluctua- 
tion relation^ for the statistics of transient current flow. 
The properties of the stochastic path integral ()2.1ip rely 
on two fundamental symmetries, namely the continuity 
equation (j2.10p and a symmetry under time reversal. The 
latter is crucial to all FRs. Our aim is thus to relate the 
functional Z, describing evolution under the influence 
of rates g = \gv,t\^ to the functional Zf, computed for 
the time- reversed rates {Tg)t = g-t, i.e., Zi, = Z\ ^ . 

Here we have defined a time reversal operator T that 
acts on "scalar" functions x = {n,g,f) as iTx)t = X-t, 
while "vectorial" functions v = (/, p, x) transform as 
iTv)t = —V-t- As shown in Appendix [Aj the action 
in Eq. (|2.1ip satisfies the symmetry 

Sg[n,p,x] = Sf^^[fn,f{p - pdnU),f{x + iPf)] 

+ P[U{nr)~U{n^r)]- (3.1) 



Next observe that n and p in Eq. p. lip are just func- 
tional integration variables. With the auxiliary relation 
for arbitrary functionals F, 



D{n,p)F[Tn,Tp] = / D{Tn,Tp) F[Tn,Tp\ 
= J D{n,p)F[n,p], 



substitution of Eq. p.ip into (|2.1ip and a shift of the 
momentum field, p -^ p + /39„t/, yields 



Z[x] = / Z?(n,p)e-*-« 



SfJn,p,T{x+i0f)] 



Pin-r)- 



We thus obtain a prototypical FR for the generating func- 
tional. 



Z[x]^Z,[f{x + tPf)] 



(3.2) 



Inserting Eq. (|3.2[) into (|2.12p . we arrive at a variant 
of the Crooks relation, first formulated by BK)^ 



Pt[fl] 



,-pn^dtj:^ui,. 



(3.3) 



where Pf, is the probability distribution computed for 
time- inverted rates Tg. Integrating over / and using the 
normalization of Pb[I], we also obtain a variant of the 
Jarzynski equality, 



(X) 



1, X = exp|/3/" dt'^fj. 



(3.4) 



Equation (|3.3p represents the most general PR relevant 
to this paper. Before turning to a discussion of its ap- 
plications, it is worthwhile to link our present classical 
formalism to the extended framework of a quantum the- 
ory of fluctuation statistics. This will be the subject of 
Sec. IIVI below. Readers primarily interested in classi- 
cal transport may skip this section and directly turn to 
Sec. El 



IV. CONNECTION TO KELDYSH APPROACH 

In this section, we show how the stochastic path inte- 
gral (|2.1ip corresponds to the /i — ;■ limit of the Keldysh 
quantum nonequilibrium functional. In this way, we will 
see how rates encoding quantum statistics, cf. Eq. (j2.5p . 
may appear as dynamical input to an effectively clas- 
sical theory of stochastic fluctuations. To be concrete, 
we focus on the example of the mesoscopic device intro- 
duced in section III Al but generalization to other setups 
is straightforward. 



A. Model 

We study a quantum dot connected to M = 2 Fermi 
liquid leads {v = 1,2) at temperatures well above the 
dot's mean level spacing and in the "open" limit, Qi, > 
1. The dimensionless conductances g,^ = 2TrhG^, with 
Gi, = R~^, describe the transparency of the contacts 
to the electrodesj^^^ We assume that the external volt- 
age Vt varies on classical time scales (such as the RC 
time of the circuit), which are large compared to quan- 
tum time scales of the problem. Technically, this means 
that terms like hdtVt can be neglected. The quantum 
nonequilibrium theory corresponds to a Keldysh func- 
tional integral)^""— 



Z = J D{<P,, cPg) exp(-5[</)„ (t)g]/h), (4.1) 

with time-dependent "classical", (j)c,t, and "quantum", 
(j)q^t, phase fields. These real- valued fields originate 
from a Hubbard-Stratonovich decoupling of the Coulomb 
interaction^^ Note that in Eq. (|4.1|) . the action S = 



Sc + 'S'tun is dimensionful, in contrast to the dimension- 
less action Sg in the classical Eq. (|2.1ip . Sc describes 
electron-electron interactions (Coulomb blockade) due to 
the charging energy Ec, 



S. 



clYc, Yq_ 



= *-^ / dt <i>qdt<i>c, 



(4.2) 



while the tunnel action contains the influence of the at- 
tached electrodes and is of Ambegaokar-Eckern-Schon 
form,— i^i^^ 



5, 



'tMn[HJc,HJq\ 



-ft E f ^' (A.e-'^Ade'^) , (4.3) 



i/=l,2 



where Kx (with X — v^A) and (j) are operators in both 
Keldysh and time space. We have cf) ~ {(t>t} with 

4>t = diag(</)+, 0,-), 0* = 0c ± 09, (4.4) 

where we employ the "contour representation" of 
Keldysh theory^ throughout. Furthermore, Ax = 
{Ax,t.t'} with a 2 X 2 matrix Ax,t,t'- For classi- 
cally varying Vt, it is convenient to pass to a Wigner 
representation)^^ whereupon matrices become functions 
of energy and time, {Ax,t,t'} -^ {Ax{t,e)}, and Tr -^ 
{2ttK)^^ j de j _ ditr, with "tr" denoting the trace in 
Keldysh space. We then have the Keldysh matri x'^^i^'^i^^ 

where nx = nx{t, e) is the electron distribution function 
in the leads {X = v) and in the dot (X = d). The 
Fermi-Dirac distribution function of the j^th reservoir is 



n^{t,e) = riF(e;K,t,7') 



1 



e(^-K)/T + 1 
where the driving voltages are 

v,.t = {~r+'vt/2. 



, (4.6) 



(4.7) 



The dot distribution function nd(i, e), however, has to be 
determined self-consistcntly. 

Noting that the product between Wigner "functions" 
under the trace in Eq. ()4.3|) has to be understood as the 
Moyal product.^ 

A*B = AB+ '-^{d.AdtB - dtAd.B) + Oift^), 
we obtain 

''^^ - 1^ 2e2^^«(l - nd.yj -(1 - 2nd,vJ ) ' 

(4.8) 
where higher-order corrections in h are neglected, 
nd,vS^) = ?^d(e — V'c), and the dynamically fluctuating 
voltage on the dot ia^ii^^ 



Vet = hdf 



(4.9) 



Substituting Eq. (|4?8)) into Eq. (|4!3]) . wc obtain 

Stun[Vc,p]^hY, fdt{l^e^P)gt{Ve) (4.10) 

with p = —2i(j)q and the rates 

5+(K) = G, j den,{e)[l-n^ySe)], 
9v{^c) = G^ / dend,y,(e)[l -rii.(e)]. 



Finally, introducing a charge-like variable through n ~ 
Vc/2Ec, adding the (n,p) representation of the charging 
term, Sc[n,p] = —hjdt pdtJi, and defining Sg[n^p\ = 
S[n^p\/h, we arrive at an action as in Eq. (|2.1ip . where 
the Hamiltonian Hg is governed by the rates 

gt{n) ^ G, f den,{e + dnU) [I - n^e)], (4.11) 



g^{n) = Gy I de[l-n^{€ + dnU)] nd(e), 
with dnU = 2Ecn. 

B. Model distributions 



This so-called "hot electron distribution" captures the 
heating of the system under the application of a voltage 
bias. We stress that Eq. (|4.12p holds only when Tee is 
short compared to the time scale for variation of Vt . (ii) 
Alternatively, one may consider a situation where the 
system is externally cooled to the ambient temperature 
T. In this case, we have 



nd(i, e) = nF(e; Voff,t,r). 



(4.14) 



(iii) In cases where the dwell time Td is smaller than the 
relaxation time Tee, the effective distribution on the dot 
is determined by the coupling to the leads rather than 
by internal relaxation. In the absence of counting fields, 
X = 0, this leads to the "double-step distribution" ob- 
tained by the weighted superposition of the two lead dis- 
tributions. 



nd{t, 



E 



9i^ 



51 +92 



-nF(e; Vi.,t,T). 



(4.15) 



However, for finite x, the situation gets more compli- 
cated in that the effective dot distribution becomes x~ 
dependent. The ensuing structures are discussed in Ap- 
pendix [B] where we also show that the FR p.2p . derived 
within a classical formalism in Sec. lIIIl stays valid in such 
a quantum-mechanical setting. 



Importantly, the rates (|4.11|) do not determine the ac- 
tion of the dot unless we specify the dot distribution func- 
tion TT-d. We here consider three different cases, all of 
which are physically relevant and conceptually interest- 
ing in their own way. Which distribution is ultimately re- 
alized depends on the ratio of two time scales, the energy 
relaxation time due to electron-electron interactions^"— 
on the dot, Tee(e) = HE^y^/e^A, and the time for escape 
into the leads, r^ = h/A{gi + 32)- Here A is the dot's 
single-particle level spacing, e is the characteristic excita- 
tion energy of particles in the system, and i?Th = S/ixh is 
the Thouless energy, where txh is the classical time scale 
before the single-particle dynamics in the dot becomes 
ergodic, e.g., the diffusion time. 

The cases considered here are as follows, (i) In the 
classical limit, h —^ 0, we observe that Tee/Td — > 0. This 
implies that strong relaxation mechanisms on the dot en- 
force an effective equilibrium Fermi distribution. 



nd{t,e) = nF(e;14ff,t,Toff,t), 



(4.12) 



where the effective chemical potential Ves and the effec- 
tive temperature Teg are determined by requiring particle 
current and energy current conservation in the dot-leads 
composite system. 
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C. Fluctuation relations 

In the remainder of this section, we address the sce- 
narios (i) and (ii) in some more detail. 

First, if the dot is cooled down to ambient temper- 
atures, the dot distribution function is given by the 
Fermi distribution (j4.14p . Substituting this function into 
Eq. (|4.11l) and doing the energy integrals, we obtain the 
rates (|2.5|) . These rates obey the detailed balance rela- 
tion (|2.4|) . which means that the externally cooled setup 
[case (ii) above] seamlessly fits into the general frame- 
work of Sees. E] and Iml Specifically, the FR dSH) and 
the derived relations in Sec. |V] below hold in full gener- 
ality. 

However, for a dot kept in isolation [case (i)], the situa- 
tion is different. Here, the temperature realized in the dot 
may differ strongly from that of the leads, which means 
that there is no uniquely specified reference temperature 
to relate to. Temperature mismatch of this type will ef- 
fectively be realized in many different circumstances: out 
of equilibrium, transport through a dissipative system 
generally leads to energy relaxation and, hence, to heat- 
ing. It stands to reason that the effective temperature 
will typically scale with the external driving parameters, 
Teff = OiVy). The resulting effective transport rates then 
no longer obey a detailed balance relation containing the 
ambient temperature as a reference scale, which in turn 
implies that the FR in Eq. (|3.3p no longer holds. 

Albeit the FR (|3.3p is violated, it is still possible to 
formulate modified FRs that contain crucial information 



about the fluctuation statistics of the system. Let us 
briefly sketch two different approaches to handling this 
situation. First, one may require that FRs categorically 
have to relate to ambient temperature, or, more gener- 
ally, to the temperatures T^, of the external leads. (In 
this part, we allow for unequal temperatures in the reser- 
voirs.) While our discussion above shows that the FR 
p.3p for particle currents is violated, it is possible to de- 
rive joint FRs for particle and energy currents that do 
apply to the heated dot yet contain the reservoir tem- 
peratures T^ only. These FRs rely on the results for 
the statistics of charge and energy transfer obtained in 
Refs. |43|68| - |7C)I . For the convenience of the reader, we 
briefly summarize the main conclusions of these refer- 
ences in the language of our paper. The main idea is to 
generalize the state space of the theory from n to (ri, e), 
where the continuous variable e represents the energy of 
the dot. It is then straightforward to derive a gener- 
alized master equation for the joint probability P{n,e), 
where the rates g^{n,e) are determined by the energy 
integrands of Eq. (|4.1ip . i.e., g'^iji) = J deg^{n,€). The 
corresponding stochastic path integral is given b y42,69/o 

Z = /"i?(n,p,e,Oe"^'''"'^'''^l, (4.16) 

Sg = - dt[pdtn + (_dte- Hg{n,p,e,^)], 

H9 = Y.j'^''{^-^^^'"'''')9t{n,e'), 



v-± 



where the time-like variable ^ is conjugate to e. In or- 
der to probe the joint statistics of particle and energy 
currentsj^i^ we couple both p and ^ to counting fields 
in the respective i^-dependent part of _ffg, 



Vt 



Pt-iXi^^t, ^t^(,t~iK 



(4.17) 



see also Eq. (|2.8p . One can then verify that the invariance 
of the action under time reversal can be effected by a 
simultaneous transformation of both counting fields, 

X. -^r(x. + i/3.K), A, ^f (A, -1/3,), 

where l3^ = T^^. In this way, the symmetry relation in 
Eq. p.l[) gets generalized and we obtain the extended FR 

Z[x., A,] = Zb[f{x, + i/3uV,), f[K - ipu)]- (4.18) 

Note that this FR involves the ambient (reservoir) tem- 
peratures T,y only. The price to be paid is that both the 
particle {xv) and the energy (Ai/) current into the utli 
reservoir have to be monitored. A quantum mechanical 
stationary version of Eq. (|4.18p has recently been dis- 
cussed in Ref . i28i. 

A second, and at this stage more heuristic, approach is 
to define time-dependent effective temperatures T* char- 
acterizing the particle exchange with the i^th reservoir 
through the logarithmic ratio of rates. 
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We now observe that the FR in Eq. (|3.2p holds provided 
one replaces the global temperature T = /3^^ by the time- 
and z^-dependent temperatures T* in Eq. (|4.19p . This 
FR contains effective (and potentially unknown) tem- 
peratures different from the ambient temperature. Fo- 
cusing on particle transport alone, it provides a reduced 
description of the nonequilibrium process. Referring to 
the applications discussed in Sec. |Vl one may employ the 
statistical information encoded in such FRs to determine 
these temperatures. This is of interest for systems where 
external biasing is expected to generate heating through 
mechanisms that are not completely understood a pri- 
ori. A precise formulation, however, is beyond the scope 
of this article and requires to carefully address several 
subtleties <ii 

In the next section, we will turn back to the general 
FR (|3.3p and discuss its applied consequences in the de- 
scription of the fluctuation statistics of nonequilibrium 
transport. 



V. APPLICATIONS 

In the present approach, the general FR p.3|) and its 
spin-offs are embedded into the formalism of the stochas- 
tic path integral. This implies a lot of freedom in explor- 
ing the role of fluctuations beyond the rigorous bounds 
implied by FRs. For example, the general Jarzynski re- 
lation p.4p states that, for a cyclic protocol, the random 
variable X averages to unity, {X) = 1. This identity 
holds under very general circumstances, and it is in this 
sense that the fluctuations of X contain more telling in- 
formation. This point has been explored in Ref. 1231 where 
we showed how the statistics of X signifles the crossover 
from near into far equilibrium situations. 

In this paper, we concentrate on the statistical infor- 
mation encoded in the FRs as such. A first aspect to 
notice is that the functional FR p.3p , which is a relation 
for the "infinitely many" variables / = {h}, contains 
information that in most applications will be excessive. 
The applied value of the identity rather lies in its poten- 
tial as a starting point for the derivation of a wealth of 
derived relations. Technically speaking, one may pass to 
these reduced identities by taking marginals in the sense 
of probability theory. Below, we discuss such reduction 
schemes on a number of examples. To keep the discussion 
concrete, we will stay with our prototypical mesoscopic 
circuit as a reference system. Generalization to other sys- 
tems should be straightforward. We stress that the FR 
p.3|) holds also in a quantum-mechanical setting, see our 
discussion above and Ref. l3ll 



A. Stationary case 

Let us begin by discussing what Eqs. p.2|) and p.3|) 
predict for the specific case of a stationary bias, Vt = V. 
More precisely, we assume a symmetric protocol Vt = V-t 



which is switched on (off) within a time Tg that is very 
short compared to the comiting time, t^ <ti t. Then, 
Vt assumes a constant value, V, during the long time 
interval |i| < t — Tg. For this voltage bias protocol, we 
have Zt, = Z. Moreover, it makes sense to consider the 
time averaged current, / = Q/2t, where Q ~ J dtit is 
the charge transmitted during the counting interval 2t. 
Taking the limit Ts/t — )■ 0, a stationary bias can then be 
described using the above formalism. 

Assuming a mesoscopic two-terminal (M = 2) setup 
for concrctcness, let us choose constant counting fields, 
Xi = — X2 = x/2, whereupon Z[x] reduces to a function 
Z{x), and differentiation w.r.t. x probes the statistics 
of /. The reduced form of Eq. I\'3.2\i then recovers the 
known FRli^sds 



Z{x) = Z{-x + i(3V), 



(5.1) 



which is equivalent to Eq. (|1.3p for the probability dis- 
tribution function P{Q). 

It is worthwhile to discuss an important consequence 
that Eq. (|5.ip entails for taking the classical limit of 
nonequilibrium quantum theories. Within the Keldysh 
approach to nonequilibrium quantum dynamics it is cus- 
tomary to associate the classical limit with a quadratic 
expansion in the quantum field, cf. Ref. (34J for a discus- 
sion of this point. However, the FR (|5.1|) implies that 
this expansion cannot be valid, unless one operates in 
a near equilibrium setting. To see this, notice that the 
counting field x couples additively to the quantum field 
(pq. Thus, a quadratic action in (j)q implies a quadratic 
S[x\- Now, consider the most general quadratic Ansatz^ 
InZ(x) = 2t[— i(/)x -I- C2X^], where we noted that dif- 
ferentiation w.r.t. X at X = yields the average current, 
(/). Consistency with Eq. (|5.1|) requires C2 = {I)/{(3V). 
For (/) r^ V, this states that the fluctuations of / (de- 
termined by the second order of the expansion in x) are 
thermal, var(/) '^ T. Thus, a quadratic expansion in 
(f>q is not capable of describing nonequilibrium noise and 
cannot be valid in general. 

More generally, Eq. (|5.ip implies constraints for the 
nonlinear coefficients describing the response of the fcth 
current cumulant {{!'')) to the voltage at order VK Since 
derivatives w.r.t. V can be traded for derivatives w.r.t. x 
in Eq. (|5.ip . different coefficients with the same order 
k+l are connected i^i^i^ The fluctuation-dissipation the- 
orem linking thermal Johnson-Nyquist noise to the lin- 
ear conductance follows from the lowest-order equation 
in this hierarchy, I + k = 2. Higher-order relations with 
l + k > 2 represent generalizations to the nonequilibrium. 
We discuss the extension of these relations to the time- 
dependent situation in Sec. IV CI 



B. Generalized Crooks relation 

We return now to a general setting with M reservoirs 
and time-dependent forces. In applications, one is often 



interested in the FCS of charge, Q[I] = J_ dtIt, trans- 
mitted at one of the terminals. (We drop the reservoir in- 
dex i> for notational simplicity.) For stationary bias, the 
probability distribution function P{Q) obeys a FR, see 
Eq. (|1.3p . and one may ask whether this relation extends 
to time-varying driving. To answer this question, we 
study the statistics of both the transmitted charge Q[I] 
and the work done on the system, Wg[/] = /_ dtltft in 
terms of their joint probability density 



PiQ,W) 



{6{Q-Q[I] 

dxqdxw 
(27r)2 ' 



)S{W-Wqm 



The second equality follows from the integral representa- 
tion of the 5- functions, where, by virtue of Eq. (|2.9p . the 
generating function Z has to be taken for the particular 
time-dependent counting field Xt = Xq ~^ Xwft- We now 
use the FR p.2p and integrate over the parameters Xq,w, 



P{Q,W) 



dXqdXv^^t{xgQ+xu,w) ^ 



'b[-Xq- iXw ~ifi)ff] 

ZbiXq+XwTf] 



(2nY 

_ f dxqdxw ,\x„Q+(x^+ii:!)w] 
J {2^Y 

^e^'^{s{Q + Q[I])5{W + Wf^[I]) 



Note that the transmitted charge and the work in the 
backward average are defined with the time-inverted 
force. In effect, we obtain a generalized Crooks relation 
coupling the transmitted charge and the work done on 
the system. 



P{Q,W) 
Pb{-Q,~W) 



aPW 



(5.2) 



which applies for time-dependent driving forces. The de- 
rived transient FR ()5.2p connects the charge to the work 
fluctuation statistics. Note that in a stationary situation, 
W (X Q, and P{Q,W) reduces to P{Q). In that case, 
Eq. (|5.2p implies the FR (|1.3p which now has the same 
physical content as Eq. (jl.ip . Turning to the generic 
time-dependent case, integrating Eq. (|5.2p over Q recov- 
ers the standard Crooks relation (|l.ip for the work dis- 
tribution function P{W). However, there is no FR for 
the reduced probability P{Q) = / dWP{Q, W) anymore, 
unless the driving force is time-independent. 



Figure [T] shows a test of the generalized Crooks relation 
(|5.2p in a numerical simulation of the master equation 
(|2.3p . Here, we have considered the asymmetric pulse 
protocol 



Vt^Vo 
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(5.3) 



where Vq is the pulse strength and 7 the pulse width. A 
relevant time scale is the RC time, J7~^ = RC/2, where 
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we choose symmetric contacts, Ri = R2 = R- We men- 
tion in passing that Pb{—Q,—W) = P{Q,—W) for an 
antisymmetric protocol hke Eq. (|5.3p . Within statisti- 
cal errors r^ the fluctuation statistics of W under a fixed 
total transferred charge Q, obtained numerically from 
4 X 10^ independent runs, agrees well with the prediction 
(El. 
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W/T 
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C. Nonlinear time-dependent transport coefficients 



FIG. 1: Numerical results for P{Q, W)/Pt{~-Q, -W) vs W/T 
for several values of the total transferred charge, taken within 
the respective window [Q — 1,Q + 1]. The straight red 
line is the result predicted by Eq. (|5.2p . Parameters: 7 — 
O.Olr, fir = 2 X 10^ Vqt = 10, and Va/T = 0.04. 



Another way to reduce the information contained in 
the transient FR for currents, Eq. p.2p . is to employ a 
series expansion of the kth current cumulant in terms of 
the driving forces fi,^t. 
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(5.4) 



Together with Eq. p.9p . this defines the time-dependent nonlinear Onsager response coefficients, 
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(5.5) 



They are symmetric under arbitrary separate permuta- 
tions of the indices {ui^ti) and [v[,t[). Due to the nor- 
malization condition Z[0] = 1, all coefficients with k ~ 
identically vanish. Moreover, for I = 0, they represent 
equilibrium correlations. Causality provides another con- 
straint: the response of the system to the driving force 
fv'.^t'. is restricted to times t > t[. This implies that 
Eq. (|5.5p must vanish whenever there is at least one 
tj > max(ti, . . . ,tk). 

Using the FR (|3.2p . we next show how for a given or- 
der k + I, different coefficients in Eq. (|5.5p are linked 
together 1^ The general relation for these time-dependent 
quantities is rather lengthy and given in App. [C] see 
Eq. (|Cip . All relations resulting for the four lowest or- 
ders (fc -I- / < 4) are also specified explicitly in App. [C] 
In experimental applications, it is often more useful to 
probe Fourier coefficients, and we provide here the spec- 
tral decomposition of the relations in App. [Cl We use the 
Fourier convention L{t) = ^ X^nez e~""*/'^L(n) (for ah 
time arguments), where the series terminates at the de- 
sired accuracy level. Since the coefficients (|5.5|) are real 
functions, we have L{~n) ~ L*{n). Rather than stating 
the general formula, we now quote the result for the few 
lowest orders, fc -t- Z = 2, 3, and 4. 



/^From Eq. (jC2p . the second-order relations are 

Im4f,)(ni,n2) = 0, (5.6) 

Rei' 



1,1)/ /\ 1^ j-{2.0) f /N 

,(ni;7ii) = --LI '{ni,n^). 
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Noting that L^^'-'^' describes the linear dissipative re- 
sponse of currents to driving, while L^^'"^ is a measure 
of equilibrium fluctuations, we recognize the second re- 
lation in Eq. (|5.6p as a spectral representation of the 
fluctuation-dissipation theorem. Moreover, the symme- 
try of Eq. (j5.5p under permutations of indices implies the 
Onsagcr-Casimir reciprocity relation 



Rei'^'^!(ni;n'i) 



Re4!:^j^(n;;ni) 
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There is no condition on Im L'^'^\ Next, from Ec^. (|C3 
we flnd the third- order relations 
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For simplicity, we here suppress the Fourier indices 
(rii, n'), which can be easily restored from the respective 
reservoir indices (z/i,i/'). The last relation in Eq. (|5.8|) 



states how the leading nonlinear response of the current 
connects to the linear order of the current noise. Finally, 
in fourth order, we get 
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The utility of such relations in the characterization of sta- 
tionary transport has been emphasized before jii^i^i^ 
Here we have applied the concept of cross-relations to 
time-dependent coefficients. Such relations provide a hi- 
erarchy of benchmark criteria for time- varying transport 
measurements or numerical simulations. 



VI. CONCLUSIONS 

In recent years, fluctuation relations have been rec- 
ognized as a potent tool in the characterization of the 
fluctuation statistics of nonequilibrium systems. Here, 
we have focused on the adaption of this concept to the 
fluctuations of transient currents, the motivation being 
that transport in response to time-varying forces repre- 
sents the perhaps most direct way of probing the nonequi- 
librium physics of complex systems. We advocated the 
stochastic path integral as an optimal tool to describe the 
emerging feedback mechanism of currents inducing noise 
which in turn affects the statistics of currents. In its 
most general form, the degree of irreversibility by which 
a transport process differs from the time-reversed pro- 
cess is characterized by the functional Crooks relation 
(|3.3[) . This relation has first been stated in the seminal 
BK paperi^ The value of the present derivation primar- 
ily lies in the linkage of the exact "sum rule" (|3.3p to the 
highly flexible formalism of the stochastic path integral. 
Indeed, we argued that the fluctuation relation is overly 
general to be of much use in concrete applications. It is 
more informative to explore fluctuations around the rig- 
orous bounds imposed by these relations, or to consider 
derived relations for the statistics of time integrated vari- 
ables, or Fourier coefficients. 

The stochastic path integral also affords an interpre- 
tation as the classical limit of the quantum theory of 
nonequilibrium fluctuations. This connection is useful 



both from a conceptual and an applied point of view. 
Conceptually, it provides a connection between the quan- 
tum and the classical theory of fluctuation relations. The 
former is usually described by means of quantum projec- 
tor techniques, a theoretical language that is not straight- 
forwardly linked to classical limits. ^From an applied 
point of view, the quantum-classical correspondence is 
of value in that it allows for an integrated description 
of processes whose stochastic rates follow from a micro- 
scopic quantum dynamics. Moreover, we have demon- 
strated that the FR stays valid in a quantum-mechanical 
regime. 

In this paper, we illustrated many of the concepts 
above on the prototypical example of a mesoscopic RC 
circuit. However, the underlying theoretical framework 
is much more general in nature, and we believe that it 
can be straightforwardly adjusted to many other different 
settings. 
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Appendix A: Symmetry of the action 

Here we provide the derivation of the symmetry re- 
lation p.ip for the action (|2.1ip of the stochastic path 
integral. Using the definition of the time-reversal opera- 
tor T in Sec. Illll and recalling the definition of K,j^^t{n) in 
Eq. (|2.4p . we have 



11 



Sf\fn,fip-(3dr.U),fix + t(3f)] = - 



dt 



p-td-tn-t + l3{dnU)-td-tn-t 



Y,(l- eT[-P-t+^X.,-t-/^«..-*(n-*)]) g±_^ 



v.± 



Next we change the time integration variable, t — > — i, and the summation variable ± — ^ =p. Noting that 
J^ dt {dnU)dtn = U{nr) — U{n^r), we obtain for the above expression 



5, 



Tg 



-/3[t/K) - U{n^r)] - 



dt 



v,± 



Ptdtnt -y^ (l 

r 






Finally, we employ the detailed balance condition (|2.4p 
and arrive at Eq. p.ip . 



We now insert Eqs. (|BT|) and ((B2|) into Eq. (PTTU]) . and 
write 5, - '^(^) - '^(^) - '^(^^ 



^tun — '-'tun 



c.(2) 
"-■tun 



S^J^. The first term reads 



Appendix B: Double-step distribution 

In this appendix)^ we discuss the structure of the tun- 
nel action (|4.3p when the effective dot distribution has 
the double-stcp form in Eq. (|4.15p . We assume that the 
system is driven and monitored on classical time scales 
such that we can neglect terms like hdtVy and hdtXy 
The driving voltages Vy were defined in Eq. ()4.7p . In the 
present quantum- mechanical context, it is crucial to in- 
clude the counting fields Xi^ from the outset. The Wigner 
representation of the A^ matrices generalizing Eq. (|4.5p 
is given by^ 



A.(t,e) = 

Ad = A/-^ 



1 - 2?!^ 2e*x..,t„^ 

9i 



91+92 

where n^it^ e) is defined in Eq. 
ensures t 
given by 



-K, 



(Bl) 



ensures the normalization condition^^^ A? = 1, and is 



The factor M{t, e) 



N 



B, 



91 +92 



v/(.9i +.92)2 +4.91.925^ 



(e'x _ i)„^(i _ n^) + ie~'^ - l)n2(l - ni). 



where x = Xi ~ X2- Note the symmetry property 
B-^^ipv = B^, which is instrumental to prove FRs. 
Keeping terms up to order fidt(l)q, instead of Eq. (|4.8p . 
the Moyal product expansion yields 



c.(i) _ " \r^ 

'-'tun ^ r, / , 



w'^ 



dt 



_2e±^(2^.+x.-x„0g±^,(^^) 



with the rates 



(B3) 



9tAyc) = 

g^AVc) = 



Gi + G2 

GyLli,! 

Gi + G2 



de [1 - n^,Ve(e)] ny'{e). 



These rates have the same functional form as the sequen- 
tial tunneling rates in Eq. (|2.5p . If (pq fluctuates on clas- 



0, we see that S[^l^ has 



sical time scales, i.e., hdt(j)q 

the same structure as the Hamiltonian Hg in Eq. ()2.8p . 

Furthermore, S^J^ coincides with S^J^, except that the 
rates g^^^, are replaced by the modified rates 

GiyLly' 



~9tAyc) = 
~9-Ayc) = 



Gi + G2 

Gi/Lri^' 

Gi + G2 



de (A/'(i,e) - 1) n^yA'^ - n^'], 
de {J\f{t,e) - 1) [1 - n^,yjn^'. 



which are now complicated functions of the dynamical 
voltage Vc, see Eq. (|4.9p , of the bias voltage V, and, 

via the normalization factor TV, of the counting field 

(2) 
X- These rates, and thus S^^^, vanish in the absence 

of counting fields, since then J\f — 1. It is easy to check 
that they still satisfy the crucial detailed balance condi- 
tion (12.41). since 



e"^Ke- 



■i(p 



(B2) 



l-2n 



i',Vc-\-hdt4>q 



2e»(x.-H20,)„ 



2e-'(x.+20,)(l _ „^_^j _(i _ 2n,,y^. 



,) 



where n^yie) = n^i^e + Vc). Notice that the energy ar- 
guments in the diagonal elements are shifted by +hdt4>q- 



r.Ay) 



^-p{V^-V^ + V^,) 



:(3) 



Finally, 5"^^;^,; = -fi(G'i + G2) / dtde {M{t, e) - 1)/2, which 
does not contribute to the field dynamics but has to be 
retained for calculating the current cumulants. 
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With the fields i/)^ in Eq. (|4.4p . we now observe that 
5tun satisfies the symmetry property 



5'tun [0t , (j^t : Xv,t 
•Jtun 



-t+tR/3/2' ^-t-ihP/2 



,-Xi^-t +il3V„,-t 



which expresses the time reversal invariance of S'tun- This 
relation, together with the invariance of Sc in Eq. 



under the replacement (j)f^ — ^ ~'f'-t±ih8/2' when inserted 
in the Keldysh generating functional (|4.ip . leads to the 
quantum generalization of the classical FR p.2p . We 
notice that in the limit S — >■ 0, the replacement 0^ — > 
— 0_ , ^„ ,„ in Eq. (|B4p is equivalently written as 

09,t ^ -(j),^-t+i(3Vc,-t/2 + Oih), 



V, 



c,t 



V,. 



o{h), 



where (/>, = ((/>+ - (f>-)/2 and V^ = ndt{(f>+ + (f>~)/2. 



Equation (|B4p thus recovers the classical relation ()3.ip 
but allows us to extend the validity of the general FR 
(B4) p.2p . which we obtained in Sec. IIIII from the classical 
generating functional, to the more general setting of the 
Keldysh quantum generating functional. 



Appendix C: Relations betw^een time-dependent 
response coefficients 

In this appendix, we provide some details concerning 
Sec. IV CI The general FR (|3.2I) for currents under time- 
dependent driving forces /^,f implies that derivatives 
w.r.t. forces can be exchanged for derivatives w.r.t. count- 
ing fields. The time-dependent Onsager coefficients ()5.5|) 
can thereby be written as 



i[^:'U;.;,...,.;(^i>--->i^;il>---.0 = n 



n h/^ 



J ' J 



Therefore we obtain a connection between coefficients with the same order k + 1, 



J ■ 3 . 



\nZ[x]\ 



x=/=o- 






(CI) 



+ /3'-iyL('+'-^'^^ _, , ,(-ti,...,-E,...,-t;;-t') 

^ / ' vi,...,Vk,i'i,...,y'^,...v[;i^'-\ J-' 'J' ' " ]' 

J = l 

I 
, al-2 V^ j(k+l-2.2) , _p T, 

^ Z-^ !^l,---,!^fc,i'J,---,P',---,P'„,---,i'i';i'',i'^i^ f-li---! f-jj-'-i f-m 

"I ^ ^!.i',...,!.fc;z.;,...,!.,'(^^l' • ■ • ' ^*fc; ^*i' ■ • ■ ' ^*i)' 



''h ''ji ''ml 



(1,0), 



where v', (resp. t'A means that i^' (resp. i') is missing in the string of indices (resp. time arguments). 

For instance, the four lowest-order relations resulting from Eq. (jCip are as follows: To first order, L\!''"'{t) = 
— lL ' (— i). The second-order result is 

4f,)(ii,i2) = 4r.l(-^i'-i2)^r4^;°l, (C2) 



vi\v{ 



I'l,!^! 



i'i:i'i 



I'l.l'i 



TL 



(1,1) 



where the time reversal operator T inverts all time arguments when acting on a function. For fc + / = 3, we find the 
relations 



TL 



(24) 



r(3,0) ^ _fr{3,0) r(2,i) ^«fr(3.") , ' 

-^1/1,1^2,1^3 1^1, 1^2, 1^3' Ui,l^2;l^'i ^1/1,1/2,1^1 1^1, 1^2; l^i' 



(C3) 



L 



(1,2) 



Finally, the fourth order produces the following relations: 



^2^,^(3,0) ,„^fM2,l) ,+i(2,l) \ 



TL 



1,2) 



r(4,0) 

1^1, 1^2, 1^3,1^4 



fr (4,0) 

Z-*! ,f2 -l-'3,?-'4 ' 



L 



(3,1) 



1/1,1/2,1/3:1/; 



_^y^(4,0) ^„y^(3,l) 

'^ Z/l,l/2,f3,^^l !^l-^'2-f3;t'i ' 



(C4) 



L 



(2,2) 



,i/2;i/i,i/2 '^ 



(4,0) 



(3,1) 



(3,1) 



, , +BT[V-'' , , +L 



TL 



(1,3) 



-(4,0) 



1^1 : f^ , i/'o , i/o '^ fi . 1/1 , Z/'A . fo '^ V fi . i/A , Z/*., ; Z/*; z/i , Z/*; . 1/4 ; Z/'A z/i . f; , z/A : i/4 



-(3,1) 



-(3,1) 



(2,2) 
1/1,1/2; 1/^,1/2' 

(3,1) 



\ l-"! iZ/i ;i'2 1^3 



(2,2) 



-(2,2) 



L' ■', , ,+L , , , 

z/i , z/2 ; z/^ , 1^3 yi . 1/3 ; 1/^ , 1^2 



Ti 



(1,3) 
'Zv*! ;Z''J ,z//, ,1/3 ■ 



13 



^ C. Bustamante, J. Liphardt and F. Ritort, Physics Today 

58, 43 (2005). 
^ G.N. Bochkov and Yu.E. Kuzovlev, Zh. Eksp. Teor. Fiz. 

76, 1071 (1979) [Sov. Phys. JETP 49, 543 (1979)]. 
3 G.N. Bochkov and Yu.E. Kuzovlev, Physica 106A, 443 

(1981). 
* R.J. Harris and G. Schiitz, J. Stat. Mech. P07020 (2007). 
^ E.M. Sevick, R. Prabhakar, S.R. Wilhams, and D.J. Sear- 

les, Annu. Rev. Phys. Chem. 59, 603 (2008). 
® U.M.B. Marconi, A. Puglisi, L. Rondoni, and A. Vulpiani, 

Phys. Rep. 461, 111 (2008). 
'' M. Esposito, U. Harbola, and S. Mukamel, Rev. Mod. 

Phys. 81, 1665 (2009). 
® L.D. Landau and E.M. Lifshitz, Statistical Physics, Part 

1, 3rd edition (Elsevier Butterworth-Heinemann, 1980). 
^ C. Jarzynski, Phys. Rev. Lett. 78, 2690 (1997). 
^° G.E. Crooks, Phys. Rev. E 60, 2721 (1999). 
" G.E. Crooks, Phys. Rev. E 61, 2361 (2000). 
^^ D. Andrieux and P. Gaspard, Phys. Rev. Lett. 100, 230404 

(2008). 
^^ M. Campisi, P. Talkner, and P. Hanggi, Phys. Rev. Lett. 

102, 210401 (2009). 
" D.M. Carberry et al, Phys. Rev. Lett. 92, 140601 (2004). 
" V. Blickle, T. Speck, L. Helden, U. Seifert, and C 

Bechinger, Phys. Rev. Lett. 96, 070603 (2006). 
^^ J.R. Gomez-Solano, A. Petrosyan, S. Ciliberto, R. 

Chetrite, and K. Gawedzki, Phys. Rev. Lett. 103, 040601 

(2009). 
" N. Gamier and S. Ciliberto, Phys. Rev. E 71, 060101 

(2005). 
^^ S. Schuler, T. Speck, C. Tietz, J. Wrachtrup, and U. 

Seifert, Phys. Rev. Lett. 94, 180602 (2005). 
^^ S. Nakamura et al, Phys. Rev. Lett. 104, 080602 (2010). 
^° Y. Utsumi et al, Phys. Rev. B 81, 125331 (2010). 
2^ R.C. Lua and A.Y. Grosberg, J. Phys. Chem. B 109, 6805 

(2005); G.E. Crooks and C. Jarzynski, Phys. Rev. E 75, 

021116 (2007). 
^^ C. Jarzynski, Phys. Rev. E 73, 046105 (2006). 
^^ A. Altland, A. De Martino, R. Egger, and B. Narozhny, 

preprint [arXiv: 1005.4662V 1 . 
■^* M. Vanevic, Yu.V. Nazarov, and W. Belzig, Phys. Rev. 

Lett. 99, 076601 (2007). 
^^ A.G. Abanov and D.A. Ivanov, Phys. Rev. Lett. 100, 

086602 (2008). 
^® J. Zhang, Y. Sherkunov, N. dAmbrumenil, and B. 

Muzykantsku, Phys. Rev. B 80, 245308 (2009). 
2^ A. Silva, Phys. Rev. Lett. 101, 120603 (2008). 
^® D. Andrieux, P. Gaspard, T. Monnai, and S. Tasaki, New 

J. Phys. 11, 043014 (2009). 
^^ P. Talkner, E. Lutz, and P. Hanggi, Phys. Rev. E 75, 

050102(R) (2007). 
^° M. Campisi, P. Talkner, and P. Hanggi, preprint 

arXiv:1003.1052vl. 
^^ M. Campisi, P. Talkner, and P. Hanggi, preprint 

larXiv:1006.1542t /l. 
^■^ A. Altland and B.D. Simons, Condensed matter field the- 
ory, 2nd edition (Cambridge University Press, Cambridge, 

2010). 
^^ Yu.V. Nazarov and Ya.M. Blanter, Quantum Transport: 

Introduction to Nanoscience (Cambridge University Press, 

Cambridge, 2009). 



A. Kamenev, in Nanophysics: Coherence and Transport, 

Les Houches session LXXXI, edited by H. Bouchiat, Y. 

Gefen, S. Gueron, G. Montambaux, and J. Dalibard (El- 
sevier, New York, 2005). 
^^ A. Kamenev and A. Levchenko, Adv. Phys. 58, 197 (2009). 
^® I.L. Aleiner, P.W. Brouwer, and L.I. Glazman, Phys. Rep. 

358, 309 (2002). 
^^ G. Schon and A.D. Zaikin, Phys. Rep. 198, 237 (1990). 
^*^ A. Ahland and R. Egger, Phys. Rev. Lett. 102, 026805 

(2009). 
^^ R. Kubo, K. Matsuo, and K. Kitahara, J. Stat. Phys. 9, 

51 (1973). 
*° K. Kitahara and H. Metiu, J. Stat. Phys. 15, 141 (1976). 
^^ P. Hanggi, Z. Phys. B 31, 407 (1978). 
"^ S. Pilgram, A.N. Jordan, E.V. Sukhorukov, and M. 

Biittiker, Phys. Rev. Lett. 90, 206801 (2003). 
^^ J. Tobiska and Yu.V. Nazarov, Phys. Rev. B 72, 235328 

(2005). 
^^ L.S. Levitov, H.-W. Lee, and G.B. Lesovik, J. Math. Phys. 

37, 4845 (1996). 
*^ H. Forster and M. Biittiker, Phys. Rev. Lett. 101, 136805 

(2008). 
'^'^ K. Saito and Y. Utsumi, Phys. Rev. B 78, 115429 (2008). 
■'^ K. Saito and Y. Utsumi, Phys. Rev. B 79, 235311 (2009). 
*^ S. Gustavsson, R. Leturcq, M. Studer, I. Shorubalko, T. 

Ihn, K. Ensslin, D.C. DriscoU, and A.G. Gossard, Surf. Sci. 

Rep. 64, 191 (2009). 
''^ T. Fujisawa, T. Hayashi, R. Tomita, and Y. Hirayama, 

Science 312, 1634 (2006). 
^^ N.G. Van Kampen, Stochastic Processes m Physics and 

Chemistry, 3rd edition (Elsevier, Amsterdam, 2007). 
^^ A. Kolomeisky and M.E. Fisher, Annu. Rev. Phys. Chem. 

58, 675 (2007). 
^^ T. Schmiedl and U. Seifert, J. Chem. Phys. 126, 044101 

(2007). 
^^ V. Mustonen and M. Liissig, PNAS 107, 4248 (2010). 
^^ R. van Zon and E.G.D. Cohen, Phys. Rev. Lett. 91, 110601 

(2003). 
^^ R. van Zon, S. Ciliberto, and E.G.D. Cohen, Phys. Rev. 

Lett. 92, 130601 (2004). 
''^ T. Taniguchi and E.G.D. Cohen, J. Stat. Phys. 130, 633 

(2007). 
^^ To keep the notation simple, we omit explicit time argu- 
ments whenever possible, rit — >■ n, etc. 
^^ Yu.V. Nazarov, Ann. Phys. (Leipzig) 16, 720 (2007). 
'^^ P.C. Martin, E.D. Siggia, and H.A. Rose, Phys. Rev. A 8, 

423 (1973). 
^° H.K. Janssen, Z. Phys. B 23, 377 (1976). 
"^ C. De Dominicis, J. Phys. (Paris), CoUoq. 37, 247 (1976). 
^^ L. Onsager and S. Machlup, Phys. Rev. 91, 1505 (1953). 
®'^ L.F. Cugliandolo, D.S. Dean, and J. Kurchan, Phys. Rev. 

Lett. 79, 2168 (1997). 
^''^ V.Y. Chernyak, M. Chertkov, and C. Jarzynski, J. Stat. 

Mech. P08001 (2006). 
®^ U. Sivan, Y. Imry, and A.G. Aronov, Europhys. Lett. 28, 

115 (1994). 
'''' Ya.M. Blanter, Phys. Rev. B 54, 12807 (1996). 
®^ B.L. Altshuler, Y. Gefen, A. Kamenev, and L.S. Levitov, 

Phys. Rev. Lett. 78, 2803 (1997). 
^^ S. Pilgram, Phys. Rev. B 69, 115315 (2004). 
^^ M. Kindermann and S. Pilgram, Phys. Rev. B 69, 155334 



14 

(2004). ''^ The slight mismatch between the slopes obtained from nu- 

T.T. Heikkila and Y.V. Nazarov, Phys. Rev. Lett. 102, merics and from the FR visible in Fig. [T] is caused by the 

130605 (2009). finite size of the window used for the numerical sampling of 

The temperatures T^ in Eq. (|4.19|l may still depend on Q, which in turn is necessary to have reasonable statistical 

the variable n. In a quasi-stationary regime with dtfi « 0, efficiency, 
we have n ~ n(V^,Ty), and a precise FR involving the 
effective temperatures T* can be stated. 



